Ionic Capillary Evaporation in Weakly Charged Nanopores 
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Using a variational field theory, we show that an electrolyte confined to a neutral cylindrical 
nanopore traversing a low dielectric membrane exhibits a first-order ionic liquid- vapor pseudo phase 
transition from an ionic-penetration "liquid" phase to an ionic-exclusion "vapor" phase, controlled 
by nanopore modified ionic correlations and dielectric repulsion. For weakly charged nanopores, 
this pseudo transition survives and may shed light on the mechanism behind the rapid switching of 
nanopore conductivity observed in experiments. 
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Electrolytes near charged surfaces are omnipresent 
in soft-matter (charged colloidal suspensions, polyelec- 
trolyte mixtures), biology (DNA, proteins and cell mem- 
branes) [T] and nanofiltration (ion channels) |3] . For 
bulk ionic fluids, the existence of an ionic first-order 
liquid-vapor (L-V) phase transition is now well estab- 
lished [3] with the driving mechanism being the competi- 
tion between the short range steric repulsion and the at- 
tractive correlations arising from the long range Coulomb 
interaction. However, for aqueous bulk electrolytes com- 
posed of conventional inorganic salts, the critical tem- 
perature is well below freezing (T c bulk w 50 K), so that 
it can be reached experimentally only with special liq- 
uids @J. When electrolytes are in contact with low di- 
electric, and possibly charged, mesoscopic bodies, inter- 
actions between mobile ions and the body surface come 
into play and strongly modify ion-ion interactions [5J [B] . 
More generally, the influence of confinement on L-V tran- 
sitions is of broad fundamental and technological inter- 
est [7]. 

In this paper, we show that when an electrolyte is con- 
fined to a neutral or weakly charged cylindrical nanopore, 
and in thermal and chemical equilibrium with an external 
salt solution reservoir, a novel type of ionic liquid-vapor 
pseudo transition occurs for conventional electrolytes at 
room temperature, in contrast to the bulk, and within 
the experimental salt concentration, pore size, and pore 
wall surface charge density range. This pseudo phase 
transition (occuring in a quasi one dimensional infinite 
system) presents parallels with capillary evaporation of 
water in hydrophobic nanopores |7j. The driving mech- 
anism is a competition between the enhanced screening 
with ionic concentration of the dielectric repulsion and 
the increase of the surface tension associated with the de- 
formation of the ionic cloud, as sketched in Fig. [I] (steric 
interactions do not seem to play an important role and 
therefore we adopt the point ion approximation). By 
using a field-theoretic variational approach, ionic corre- 
lations and polarization charge effects can be taken into 
account non-perturbatively (Sj [9] . In the bulk, if l/<Zj 
is introduced as a cutoff in momentum space, this ap- 
proach correctly predicts the existence of a first-order 




FIG. 1: Sketch of a cylindrical nanopore (radius a) filled with 
counterions and coions. The surface charge density is a s , k v 
and Kb are screening parameters in the pore and in the bulk 
(cylindrical coordinates r,z). 



ionic L-V transition |10j . In neutral nanopores, we also 
find a (pseudo) transition between a high concentration, 
conducting ionic liquid phase and a very low concentra- 
tion insulating ionic-exclusion vapor phase; for weakly 
charged nanopores, however, the pseudo transition is to 
a low conductivity counterion-only vapor phase, where 
coions are nearly entirely excluded from the pore and 
due to electroneutrality the counterion concentration is 
fixed almost entirely by the surface charge density. We fi- 
nally propose that the underlying mechanism controlling 
conductivity fluctuations as a function of pH and diva- 
lent ion concentration in certain artificial and biological 
nanopores [TT1H5] is a manifestation of the pseudo tran- 
sition proposed here in a weakly charged nanopore. 

We consider an electrolyte at T — 300 K of dielectric 
permittivity e w = 78, confined in a cylindrical nanopore 
of radius a, length L, and surface charge density cr s ; 
the space outside the pore is salt free with a dielec- 
tric permittivity e m < e w (Fig. [I]). The electrolyte is 
in contact with an external ion bulk reservoir at the 
end boundaries of the pore, which fixes the fugacity of 
ions inside the pore according to chemical equilibrium, 
Aj = e Mi /A 3 = Aj,!,, where fa is the chemical potential of 
ion i = 1 . . . a (energies are in units of the thermal energy 
ksT = 1//3) and A is the De Broglie wavelength of an 
ion. Although ions interact through the bare Coulomb 
potential, Vq(t, r') = £b/\ t — r '| in a bulk electrolyte, 
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where £b = /3e 2 /(47re w ) w 0.7 nm, dielectric jumps at the 
nanopore boundaries yield a modified Coulomb potential 
uc(r, r') obeying — Ve(r)Vvc(r, r') = /3e 2 <5(r — r'). After 
performing a Hubbard-Stratonovich transformation and 
introducing a fluctuating field 4>(r), the grand-canonical 
partition function for ions in the nanopore can be written 
as Q = fV<j> e~ H ^/Q c where the Hamiltonian is [HI E] 
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qi is the ion valency, Qc = — ^trln^c, c(r) = cr s <5(r — a), 
and Aj = Aie 9 ^^ )/ 2 . The electrostatic potential is given 
by i(0(r)). 

Evaluating the partition function using Eq. ([I]) is in- 
tractable due to non-linear terms, and we use the Gaus- 
sian variational method, which consists in computing 
the variational grand-potential Q v — + (H — Hq)q, 
where the expectation value is evaluated with a varia- 
tional Gaussian Hamiltonian 
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and f^o = — |trln(w /uc)- We then extremize Q v with 
respect to the variational functions, namely the Green's 
function vo(r,r') and the electrostatic potential 0o(r). 
This yields two intractable coupled non-linear differen- 
tial equations, similar in form to a generalized Poisson- 
Boltzmann equation for </>o( r ) and a generalized Debye- 
Hiickel (DH) one for vo(r,r') [5J |S]. To make progress 
we consider the restricted case of a constant cj>o in the 
pore [2] and take vo(r,r') as the solution of [— Ve(r)V + 
e(r)K 2 (r)]«o(r, r') = f3e 2 5(r — r') with a screening param- 
eter, /c(r) = k v 0(a — r) [2J[B]. We are then left with two 
variational parameters, the effective Donnan potential, 
00) and the DH parameter in the pore k v [5]. 

The variational grand-potential becomes f2„ = 
-pira 2 L + j s v 2iraL, where p = J2i \e q ^ BR "l 2 - k 3 J(24it) 
is the pressure of a bulk electrolyte with screening pa- 
rameter k v |14j and 7^ a surface contribution given by 
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where (. . .) is the average over the nanopore and 
v (r, r') = f B e- K -l r '- r l /|r' - r| + $u (r, r') with 

<Su (r,r; «„) = / dfc V" F m (/s; K„)I 2 a (xr) (4) 

1/0 m>0 

where >r 2 = /c 2 + k 2 , ^2 means that the term to = is 
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and I m and K m are modified Bessel functions [T51 [TB] . 
The first term in Eq. ([3| is the electrostatic energy of the 
surface charge, the second term is a depletion term and 
the last term is the cost of ionic cloud deformation. 

In the following, we consider low to moderate ion con- 
centrations in the bulk reservoir such that the physi- 
cal minimum of Q v in the bulk (for a — > 00) is given 
by the DH result E], X lyb = p ij6 e -<?,W B /2 where 
k 2 = A-kIb J2t=i QiPi-b- The ion concentration aver- 
aged over the pore section is (pi) = Pi : b^ie~ qi ^° with 

= (e~ q i w ( r ^ 2 ) where the potential w(r) incorporates 
the solvation energy due to the ionic cloud and polariza- 
tion charge-ion interactions 



w(r) = (K b - k v )£ b + 8v (r,r; k v ) 



(6) 



and $i(r) = q 2 w(r)/2 + is the potential of mean 
force of ion i in the nanopore [9] . 

For a symmetric electrolyte (q± = q, T± = T), by 
extremizing il v with respect to (f>o, one obtains the elec- 
troneutrality condition, o~ s = qpbaT sinh (qcfio), which, in 
the mean-field Poisson-Boltzmann limit (w(r) = 0) leads 
to the usual Donnan potential [2J. By injecting the so- 
lution for in Q v , we are left with a single variational 
parameter n v (pb, a, e m Aw)- The averaged partition coef- 
ficients of countcrions and coions arc 



Pb 



(7) 



In a neutral nanopore, <p — and k± = T. For mono- 
valent ions (q = 1) and low enough membrane permittiv- 
ities (e m < 5), pore radii (a < 1 nm) and bulk ion concen- 
trations, the variational grand-potential Q v (k v ) exhibits 
a minimum at K^a ~ 5 x 10~ 2 . By increasing the reser- 
voir concentration this minimum jumps discontinuously 
to a finite value k% with < < Kj,, and the pore 
undergoes a first-order pseudo phase transition from an 
ionic- exclusion state to an ionic-penetration one. As il- 
lustrated in Figs. [2^, and[2j} for e m = 2, where k is plot- 
ted, respectively, vs pb and a, for a < a* — 0.989 nm, a 
critical pore radius, a jump occurs for a specific coexis- 
tence value of pb larger than pi = 0.074 mol/L. At the 
transition, there is phase coexistence, Q, v (k%) = Q v (k%), 
where k^/k^ ~ 10. When increasing pb beyond this 
coexistence value, the ionic-exclusion stable state be- 
comes metastable and the pore becomes penetrable to 
ions (Fig. [2^,) . For a = a* , the pseudo transition is con- 
tinuous. 

Because of the strong ion depletion near the pore sur- 
face, an instructive analogy with capillary evaporation of 
water in hydrophobic nanopores favoring the vapor phase 
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FIG. 2: Partition coefficient k = {p)/pb inside a neutral nanopore (q = 1, = 78) vs (a) the bulk concentration p ft 

for three pore radii with metastable branches (thin solid lines) and window (vertical lines) for a = 0.84 nm; and (b) the pore 
radius a for, from left to right, pt = 0.7,0.3,0.156,0.08 mol/L. Dotted (grey/red) lines show metastable (unstable) branches, 
light grey/blue lines (guide for the eye) are the "boiling point" curve (bottom) and the "dew point" curve (top) and the dot is 
the critical point, pi — 0.074 mol/L, a* = 0.989 nm. (c) Phase diagram for various e m (e w = 78). Critical lines correspond to 
phase boundaries between the ionic-penetration phase (L, above) and the ionic-exclusion one (V, below). 



can be drawn [7]. For ionic fluids, in contrast to water, 
we find that the critical temperature increases with de- 
creasing pore radius. Whereas the critical temperature 
T c of the L-V transition in bulk electrolytes is extremely 
low (reduced temperature T bulk = 2a,i/£B — 0.05 with 
dj ^ 0.1 nm the ion radius), the phase diagram in Fig.[2j; 
shows that at room temperature, the exclusion pseudo 
phase transition takes place for experimentally accessi- 
ble parameter values (e w = 78, e m = 1 to 4, and re- 
duced temperature TP orc = 2a/ £b — 1): the coexistence 
lines separate the L-, above the line, from V-state, be- 
low the line. Beyond the critical point (a*(e m ),/3^(e m )), 
we enter the "fluid" phase where the pseudo transition 
disappears and is replaced by a smooth crossover (as in 
slit pores (2j [9j [XT] ). Using a simplified self-consistent 
approach, Dresner also found a first-order pseudo phase 
transition for electrolytes confined in neutral spherical 
pores [17] . but Yaroschuk may have wrongly argued 
that this is an artefact arising from the use of the self- 
consistent DH equation [2] . When increasing e m and thus 
decreasing the repulsive surface polarization charge, the 
parameter range where the phase separation is observable 
is considerably reduced. Hence, for neutral pores the ion 
penetration pseudotransition is driven by the competi- 
tion between the last two terms of Eq. ([3|, which favor 
respectively high and low k v . 

What happens for slightly charged nanopores? In 
Fig. [3^i are plotted the partition coefficients of counte- 
rions k+ and coions fc_ and the variational Donnan po- 
tential </>o vs for e m = 2, a — 0.84 nm and a weak 
surface charge, <r s — —8.6 x 10 -5 C/m 2 . One observes 
that the discontinuous pseudo phase transition survives 
and phase coexistence occurs at pb — 0.17 mol/L, a value 
slightly lower than for the neutral case (0.18 mol/L). At 
this coexistence value, 4>o a l so exhibits a jump, and for 
larger pb, (f>Q ~ and we recover the neutral case. A slight 
difference between k+ and fc_ remains due to electroneu- 



trality, rewritten using Eq. ([T]) as k + — /c_ = 2|cr s |/ (qp^a). 
For smaller pb, k- <C k + while k + rapidly increases. 
We thus reach a low surface charge density counterion- 
only regime, with k + ~ 2\a s \/qpba (Fig. [3^,, dotted 
line) and /c_ = T 2 /k+ [9]. In this regime, ~ 
2\a s \/(qa), is independent of pb and determined solely by 
global electroneutrality, and coions are excluded mainly 
by dielectric repulsion and not charge. Hence the va- 
por phase is no longer an ionic-exclusion phase but a 
weak ionic-penetration one, governed by the weak sur- 
face charge density. For larger surface charge densities, 
|cr s | > 10~ 3 C/m 2 , the transition disappears. The phase 
diagram for several values of a s is illustrated in Fig. [3]d. 
Increasing the surface charge favors the L phase by re- 
ducing the coexistence line and shifting the critical point 
towards smaller a and larger pb (because o~ s increases 7^ , 
but due to screening, this enhancement of 7^, decreases 
with K v ). Comparison of Figs. [2J: and[3}D clearly shows 
that an increase in |er s | plays qualitatively the same role 
as an increase in e m , i.e. a decrease in the repulsive po- 
larization charge density. 

Can the trace of this first-order pseudo transition be 
observed experimentally in finite sized open-ended pores? 
By studying ionic conductivity in nanopores (a ~ 1 nm, 
L ~ 5 /im) produced in polyethylene terephthalate mem- 
branes [TT], Lev et al. observed three regimes: a high 
conductivity (HC) regime, a low conductivity (LC) one, 
and a two-state high/low conductivity (HC/LC) regime 
where the conductivity switches rapidly between both 
states. Interestingly, the dynamic characterization of the 
HC/LC regime, performed by measuring the ratio of resi- 
dent times, leads to the identification of three parts: HC 
more stable, HC and LC at coexistence, and LC more 
stable. We argue that one possible mechanism for these 
current fluctuations on the scale of seconds is that the 
system is close to the phase coexistence presented here, 
where the HC and LC regimes are identified, respectively, 
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FIG. 3: (a) Partition coefficients of coions fc_ and counterions k+ (left) and effective Donnan potential 4>o (right) vs. pj, for 
charged pores (e m = 2, a = 0.84 nm, and fixed surface charge density a s = —0.9 x 10~ 4 C/m 2 ). The dotted line is the 
counterion-only approximation for counterions. (b) Phase diagram (similar to Fig. j2ji) for several values of \<r s \. Inset: k± vs. 
\a s \ (pb = 1 mol/L, a = 0.617 nm) showing the window (horizontal arrow) with stable (solid lines) and metastable branches 
(dashed lines) ; the predicted HC to LC conductivity ratio oc ^ kf / kY — 5 is close to that seen in experiments [TT] (a value 
difficult to explain via water capillary evaporation [7]). (c) Logarithm of the ratio of resident times in L and V states of the 
inset of (b), tl/tv = exp[f2(fcjf) — fi(ft^)] vs. \a s \ (pore length L — 5 fim [TT]). 



with the ionic-penetration "liquid" and ionic-weak pene- 
tration "vapor" phases, and the HC/LC regime with the 
window of <j s for which metastable branches exist (inset 
of Fig. [3}d). For finite open-ended nanopores, we expect 
the L-V pseudo transition to be rounded, with phase sep- 
aration and hysteresis replaced by two-state fluctuations 
between pseudo-stable or metastable phases [7] . For sur- 
faces carrying carboxylic acid groups, \a s \ is an increas- 
ing function of pR [TS] (for more details see [H]), and 
Lev et al. jTTJ showed that the HC/LC regime only ex- 
ists within a narrow pH window. At high pK, |<r g | is 
high and only the HC state exists; at low pH, \a s \ is low 
and only the LC state exists. Within this "fluctuation" 
pH window, using the two-state approximation [7J, our 
model yields trends for ratio of resident times vs. |ct s | 
(Fig. [3J;) in qualitative agreement with the experimen- 
tal ones vs. pH (Fig. 4 of [UJ). Similarly, increas- 
ing the concentration of trace divalent cations decreases 
the bare negative surface charge [Tj5] and the nanoporc 
can likewise switch from the HC to the LC state [TT] . 
Observations revealing a strong correlation between a s 
and conductivity fluctuations [12] can also be explained 
via the charge-regulation mechanism [TH] , because a s de- 
pends in a self-consistent way on local solution charac- 
teristics |16j . Although other mechanisms have been pro- 
posed to explain nanopore conductivity fluctuations (gas 
bubbles [13] or salt occlusions [20]), we believe that the 
one proposed here provides a natural explanation for a 
host of experimental trends [TJJ [H] and deserves further 
detailed investigation [16) . One open question concerns 
the role of ion-pairing; such effects appears naturally in 
a 2nd order variational approach and are currently under 
study. bUsing the Fisher and Levin [3] bulk approach as 
a guideline, we expect higher order corrections to lead to 
quantitative, but not qualitative, changes to our results. 



To corroborate our predictions it would also be interest- 
ing to perform simulations, although it may be extremely 
difficult to properly include dielectric discontinuities and 
reach sufficiently long time scales (~ 1 s). 
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